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PULSE OXIMETER WITH SEPARATE ENSEMBLE AVERAGING FOR 
OXYGEN SATURATION AND HEART RATE 

CROSS-REFERENCES TO RELATED APPLICATIONS 
5 [0001] NOT APPLICABLE 

STATEMENT AS TO RIGHTS TO INVENTIONS MADE UNDER 
FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT 
[0002] NOT APPLICABLE 

10 

REFERENCE TO A "SEQUENCE LISTING," A TABLE, OR A COMPUTER 
PROGRAM LISTING APPENDIX SUBMITTED ON A COMPACT DISK 
[0003] NOT APPLICABLE 

1 5 BACKGROUND OF THE INVENTION 

[0004] The present invention relates to oximeters, and in particular to ensemble averaging 
of pulses in a detected waveform from a pulse oximeter. 

[0005] Pulse oximetry is typically used to measure various blood chemistry characteristics 
including, but not limited to, the blood-oxygen saturation of hemoglobin in arterial blood, the 

20 volume of individual blood pulsations supplying the tissue, and the rate of blood pulsations 
corresponding to each heartbeat of a patient. Measurement of these characteristics has been 
accomplished by use of a non-invasive sensor which scatters light through a portion of the 
patient's tissue where blood perfuses the tissue, and photoelectrically senses the absorption of 
light in such tissue. The amount of light absorbed at various wavelengths is then used to 

25 calculate the amount of blood constituent being measured. 

[0006] The light scattered through the tissue is selected to be of one or more wavelengths 
that are absorbed by the blood in an amount representative of the amount of the blood 
constituent present in the blood. The amount of transmitted light scattered through the tissue 
will vary in accordance with the changing amount of blood constituent in the tissue and the 
30 related light absorption. For measuring blood oxygen level, such sensors have typically been 



provided with a light source that is adapted to generate light of at least two different 
wavelengths, and with photodetectors sensitive to both of those wavelengths, in accordance 
with known techniques for measuring blood oxygen saturation. 

[0007] Known non-invasive sensors include devices that are secured to a portion of the 
5 body, such as a finger, an ear or the scalp. In animals and humans, the tissue of these body 
portions is perfused with blood and the tissue surface is readily accessible to the sensor. 

[0008] N-100 . The N-100 technology, dating to around 1985, accepted or rejected pulses 
based on pulse history of the size of pulses, pulse shape, expected time to occur (frequency) 
and ratio of R/IR. 

10 [0009] In particular, the N-100 found pulses by looking for a signal maximum, followed by 
a point of maximum negative slope, then a minimum. The processing was done in a state 
machine referred to as "munch." Each maximum was not qualified until the signal passed 
below a noise threshold, referred to as a noise gate. This acted as an adaptive filter since the 
noise gate level was set by feedback from a subsequent processing step to adapt to different 

15 expected signal amplitudes. The pulses are then accepted or rejected in a "Level3" process 
which was a filter which adapts to changing signals by comparing the amplitude, period and 
ratio-of-ratios (ratio of Red to IR, with Red and IR being expressed as a ratio of AC to DC) of 
a new pulse to the mean of values in a history buffer, then determining if the difference is 
within a confidence level. If the new pulse was accepted, the history buffer was updated with 

20 the values for the new pulse. The leveB process acted as an adaptive bandpass filter with 
center- frequency and bandwidth (confidence limits) being adapted by feedback from the 
output of the filter. 

[0010] N-200 . The N-200 improved on the N-100 since it could be synchronized with an 
ECG, and included ECG filtering. The N-200 also added interpolation to compensate for 
25 baseline shift between the time of measuring the pulse maximum and minimum. The N-200 
included other filtering features as well, such as a "boxcar" filter which computed the mean 
of a varying number of signal samples. 

[0011] The N-200, after various filtering and scaling steps, applies the digitized signals to a 
"boxcar" filter, which computes the mean of N samples, where N is set by feedback from a 
30 subsequent processing step according to the filtered heart rate. New samples are averaged 

into the boxcar filter, while the oldest samples are dropped. The boxcar length (N) is used to 
set three parameters: a pulse threshold, absolute minimum pulse and small pulse. An 
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ensemble-averaging (a.k.a "slider") filter then produces a weighted average of the new 
samples and the previous ensemble-averaged sample from one pulse-period earlier. The 
samples are then passed to a "munch" state machine and a noise gate, like the N-100. An 
interpolation feature is added to the N-100 process, to compensate for changes in the baseline 
5 level. Since the minimum and maximum occur at different times, a changing baseline may 
increase or decrease the minimum and not the maximum, or vice-versa. 

[0012] "Ensemble averaging" is an integral part of C-Lock, which is NELLCOR's 
trademark for the process of averaging samples from multiple pulses together to form a 
composite pulse. This process is also known as "cardiac-gated averaging." It requires a 
10 "trigger" event to mark the start of each pulse. 

[0013] Conlon US Patent No. 4,690,126 discloses ensemble averaging where different 
weights are assigned to different pulses and a composite, averaged pulse waveform is used to 
calculate oxygen saturation. The N-100 described above is described in US Patent No. 
4,802,486. Aspects of the N-200 are described in US Patent No. 4,91 1,167 (Corenman) and 
15 No. 5,078,136 (Stone). 

BRIEF SUMMARY OF THE INVENTION 
[0014] The present invention is directed to the use of two separate ensemble averagers for 
processing a detected waveform for use in calculating oxygen saturation and a pulse rate. 
20 The ensemble averager used for calculating oxygen saturation operates on a signal which has 
been normalized, while the ensemble averager for the pulse rate calculation operates on a 
signal which has not been normalized. Note that the waveforms corresponding to both 
wavelengths must be normalized by the same quantity, such as the IR pulse amplitude, so as 
to preserve the ratio-of-ratios for oxygen saturation computation. 

25 [0015] The use of a signal without normalization for the pulse rate improves the software's 
ability to disqualify artifacts that are substantially larger than physiological pulses, such as 
motion artifact. The use of a signal without normalization for the pulse rate avoids a pulse 
being missed due to normalization. 

[0016] The metrics chosen for the two paths through the two ensemble averagers can be 
30 varied to optimize the ensemble averaging for oxygen saturation or pulse rate calculations. 
For example, a lower threshold is used for a metric to detect arrhythmic pulses when used to 
calculate pulse rate, as compared to calculating oxygen saturation. Also, a metric for a short 
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term pulse amplitude ratio will be small when motion artifact has just subsided, and this is 
given more weight in the pulse rate calculation than in the oxygen saturation calculation (the 
short-term pulse amplitude ratio is current pulse amplitude / previous pulse amplitude). 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

[0017] Fig. 1 is a block diagram of an oximetry system incorporating an embodiment of the 
invention. 

[0018] Fig. 2 is a diagram of the software processing blocks of an oximeter including an 
embodiment of the present invention. 

10 [0019] Fig. 3 is a diagram showing the creation of a composite pulse. 

[0020] Fig. 4 is a chart of the ensemble averaging performance. 

[0021] Figs. 5-7 are diagrams of state machines for updating certain variables in a 
composite pulse buffer. 

1 5 DETAILED DESCRIPTION OF THE INVENTION 

[0022] Fig. 1 illustrates an embodiment of an oximetry system incorporating the present 
invention. A sensor 1 0 includes red and infrared LEDs and a photodetector. These are 
connected by a cable 12 to a board 14. LED drive current is provided by an LED drive 
interface 16. The received photocurrent from the sensor is provided to an I-V interface 18. 

20 The IR and red voltages are then provided to a sigma-delta interface 20 incorporating the 

present invention. The output of sigma-delta interface 20 is provided to a microcontroller 22 
which includes a 10-bit A/D converter. Controller 22 includes flash memory for a program, 
and EEPROM memory for data. The processor also includes a controller chip 24 connected 
to a flash memory 26. Finally, a clock 28 is used and an interface 30 to a digital calibration 

25 in the sensor 10 is provided. A separate host 32 receives the processed information, as well 
as receiving an analog signal on a line 34 for providing an analog display. 

[0023] Design Summary The design of the present invention is intended to deal with 
unwanted noise. Signal metrics are measured and used to determine filter weighting. Signal 
metrics are things that indicate if a pulse is likely a plethysmograph or noise, such as 
30 frequency (is it in the range of a human heart rate), shape (is it shaped like a heart pulse), rise 
time, etc. A similar technique was used in the Nellcor N200, described in the background of 
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this application. The new design adds a number of different features and variations, such as 
the use of two ensemble averagers as claimed in the present invention. 

[0024] Details of the architecture are shown in the diagram of Fig. 2. This design 
calculates both the oxygen saturation, and the pulse rate, which are described separately 
5 below. 

[0025] I. Oxygen Saturation Calculation. 

[0026] A. Signal Conditioning - The digitized red and IR signals are received and are 
conditioned in this block by (1) taking the 1st derivative to get rid of baseline shift, (2) low 
pass filtering with fixed coefficients, and (3) dividing by a DC value to preserve the ratio. 
10 The function of the Signal Conditioning subsystem is to emphasize the higher frequencies 
that occur in the human plethysmograph and to attenuate low frequencies in which motion 
artifact is usually concentrated. The Signal Conditioning subsystem selects its filter 
coefficients (wide or narrow band) based on hardware characteristics identified during 
initialization. 

15 [0027] Inputs - digitized red and IR signals 

[0028] Outputs - Pre-processed red and IR signals 

[0029] B. Pulse Identification and Qualification - The low pass filtered and digitized red 
and IR signals are provided to this block to identify pulses, and qualify them as likely arterial 
pulses. This is done using a pre-trained neural network, and is primarily done on the IR 
20 signal. The pulse is identified by examining its amplitude, shape and frequency, just as was 
done in the Nellcor N-100. An input to this block is the average pulse period from block D. 
This function is similar to the N-100, which changed the upfront qualification using the pulse 
rate. The output indicates the degree of arrhythmia and individual pulse quality. 

[0030] Inputs - (1) Pre-processed red and IR signals, (2) Ave. pulse period, (3) 

25 Lowpass Waveforms from the low pass filter. 

[0031] Outputs - (1) Degree of arrhythmia, (2) pulse amplitude variations, (3) 

individual pulse quality, (4) Pulse beep notification, (5) qualified pulse periods and age. 

[0032] C. Compute Signal Quality Metrics - This block determines the pulse shape 
(derivative skew), period variability, pulse amplitude and variability, Ratio of Ratios 
30 variability, and frequency content relative to pulse rate. 
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[0033] Inputs - (1) raw digitized red and IR signals, (2) degree of arrhythmia, 

individual pulse quality, pulse amplitude variation (3) pre-processed red and IR signals, 
(4) average pulse period. 

[0034] Outputs - (1) Lowpass and ensemble averaging filter weights, (2) metrics for 

5 sensor off detector, (3) Normalized Pre-processed waveforms, (4) percent modulation. 

[0035] D. Average Pulse Periods. This block calculates the average pulse period from the 
pulses received. 

[0036] Inputs - Qualified pulse periods and age. 

[0037] Outputs - Average pulse period. 

10 [0038] El. Lowpass Filter and Ensemble Averaging - Block El low pass filters and 

ensemble averages the signal conditioned by block A, and normalized by block C, for the 
pulse rate identification. The weights for the low pass filter are determined by the Signal 
Metrics block C. The signal is also ensemble averaged (this attenuates frequencies other than 
those of interest near the pulse rate and its harmonics), with the ensemble averaging filter 

1 5 weights also determined by Signal Metrics block C. Less weight is assigned if the signal is 
flagged as degraded. More weight is assigned if the signal is flagged as arrhythmic because 
ensemble-averaging is not appropriate during arrhythmia. Red and IR are processed 
separately, but with the same filtering weights. The filtering is delayed approximately one 
second to allow the signal metrics to be calculated first. 

20 [0039] The filters use continuously variable weights. If samples are not to be ensemble- 
averaged, then the weighting for the previous filtered samples is set to zero in the weighted 
average, and the new samples are still processed through the code. This block tracks the age 
of the signal - the accumulated amount of filtering (sum of response times and delays in 
processing). Too old a result will be flagged (if good pulses haven't been detected for 

25 awhile). 

[0040] Inputs - (1) normalized pre-processed red and IR signals, (2) average pulse 

period, (3) low pass filter weights and ensemble averaging filter weights, (4) ECG triggers, if 
available, (5) IR fundamental, for zero-crossing triggers. 

[0041] Outputs - (1) filtered red and IR signals, (2) age. 
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[0042] F. Estimate Filtered Waveform Correlation and Calculate Averaging Weight - this 
uses a noise metric similar to that used in the N100 and N200 described above, and doesn't 
use feedback. The variable weighting for the filter is controlled by the ratio-of-ratios 
variance. The effect of this variable-weight filtering is that the ratio-of-ratios changes slowly 
5 as artifact increases and changes quickly as artifact decreases. The subsystem has two 

response modes. Filtering in the Fast Mode targets an age metric of 3 seconds. The target age 
is 5 seconds in Normal Mode. In Fast Mode, the minimum weighting of the current value is 
clipped at a higher level. In other words, a low weight is assigned to the newest ratio-of- 
ratios calculation if there is noise present, and a high weight if no noise is present. 

10 [0043] Inputs - (1) filtered red and IR signals and age, (2) calibration coefficients, (3) 

response mode (user speed settings). 

[0044] Outputs - averaging weight for ratio-of-ratios calculation. 

[0045] H. Calculate Saturation - Saturation is calculated using an algorithm with the 
calibration coefficients and averaged ratio of ratios. 

15 [0046] Inputs - (1) Averaged Ratio-of-Ratios, (2) calibration coefficients. 

[0047] Outputs - Saturation. 

[0048] II. Pulse Rate Calculation. 

[0049] E2. Lowpass Filter and Ensemble Averaging - Block E2 low pass filters and 
ensemble averages the signal conditioned by block A, for the pulse rate identification. The 

20 weights for the low pass filter are determined by the Signal Metrics block C. The signal is 

also ensemble averaged (this attenuates frequencies other than those of interest near the pulse 
rate and its harmonics), with the ensemble averaging filter weights also determined by Signal 
Metrics block C. Less weight is assigned if the signal is flagged as degraded. More weight is 
assigned if the signal is flagged as arrhythmic since filtering is not appropriate during 

25 arrhythmia. Red and IR are processed separately. The process of this block is delayed 
approximately one second to allow the signal metrics to be calculated first. 

[0050] The filters use continuously variable weights. If samples are not to be ensemble- 
averaged, then the weighting for the previous filtered samples is set to zero in the weighted 
average, and the new samples are still processed through the code. This block tracks the age 
30 of the signal - the accumulated amount of filtering (sum of response times and delays in 
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processing). Too old a result will be flagged (if good pulses haven't been detected for 
awhile). 

[0051] Inputs - (1) pre-processed red and IR signals, (2) average pulse period, (3) 

Lowpass filter weights and ensemble averaging filter weights, (4) ECG triggers, if available, 
5 (5) IR fundamental, for zero-crossing triggers. 

[0052J Outputs - (1) filtered red and IR signals, (2) age. 

[0053] I. Filtered Pulse Identification and Qualification - This block identifies and 
qualifies pulse periods from the filtered waveforms, and its results are used only when a pulse 
is disqualified by block B. 

10 [0054] Inputs - (1) filtered red and IR signals and age, (2) average pulse period, (3) 

hardware ED or noise floor, (4) kind of sensor. 

[0055] Outputs - qualified pulse periods and age. 

[0056] J. Average Pulse Periods and Calculate Pulse Rate - This block calculates the pulse 
rate and average pulse period. 

1 5 [0057] Inputs - Qualified pulse periods and age 

[0058] Outputs - (1) average pulse period, (2) pulse rate. 

[0059] III. Venous Pulsation 

[0060] K. Detect Venous Pulsation - Block K receives as inputs the pre-processed red and 
IR signal and age from Block A, and pulse rate and provides an indication of venous 
20 pulsation as an output. This subsystem produces an IR fundamental waveform in the time 
domain using a single-tooth comb filter which is output to the Ensemble Averaging filters. 

[0061] Inputs - (1) filtered red and IR signals and age, (2) pulse rate. 

[0062] Outputs - Venous Pulsation Indication, IR fundamental 

[0063] IV. Sensor Off 

25 [0064] L. Detect Sensor-Off and Loss of Pulse Amplitude - The Pulse Lost and Sensor Off 
Detection subsystem uses a pre-trained neural net to determine whether the sensor is off the 
patient. The inputs to the neural net are metrics that quantify several aspects of the behavior 
of the IR and Red values over the last several seconds. Samples are ignored by many of the 



algorithm's subsystems while the Signal State is not either Pulse Present or Sensor Maybe 
Off. The values of the Signal State variable are: "Pulse Present, Disconnect, Pulse Lost, 
Sensor Maybe Off, and Sensor Off. 11 

[0065] Inputs - (1) metrics, (2) front-end servo settings and ID 

5 [0066] Outputs - Signal state including sensor-off indication 

ENSEMBLE AVERAGING SUBSYSTEM 

[0067] Ensemble Averaging subsystem . The function of the Ensemble Averaging 
subsystem is to filter its input streams with a variable weighting, and output waveforms that 
are less distorted by noise or motion artifact. Reducing the degree of artifact in the filtered 
10 waveforms enables a more robust saturation or rate estimate during motion or noise. 

[0068] The Ensemble Averaging subsystem requires IR and Red inputs every sample that 
are zero-mean over the span of several pulses. 

[0069] Input samples are first IIR lowpass filtered with a weight received from the Signal 
Metrics subsystem (LPF Weight) and then stored in a one-second-delay buffer. 

1 5 [0070] It averages the zth IR and Red samples of the current one-second delayed input 
pulse with the ith samples of the previous composite pulse to form the ith samples of the 
current composite pulse. The trigger to start the beginning of a pulse period is derived from 
(in order of priority) the R WaveJDccurred input and the average period input 
(Optical _Period). The weight given to the current sample versus the corresponding sample 

20 of the previous pulse is determined by the value of Ensemble _Aver aging JVeight value 
received from the Signal Metrics subsystem. 

[0071] This composite pulse is less distorted by noise or motion artifact than the filter's 
input pulses. Figure 3 is a conceptual illustration of how pulses are averaged together to form 
the composite pulse. 

25 [0072] The subsystem receives and qualifies triggers, which should be synchronous with 
the heartbeat. The triggers are qualified R-Wave triggers from the R-wave Qualification 
subsystem when available. When R-Wave triggers are unavailable, the triggers are generated 
internally from the average period input from the Pulse Rate Calculation subsystem 
{Optical _Period). A "pulse" is considered to start at each qualified trigger and end at the 
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next qualified trigger. In this way, consecutive triggers are used to define the ensemble 
averaging period. 

[0073] Figure 4 shows a representation of the subsystem's response to a series of pulses 
corrupted by motion artifact. The vertical lines are R-Wave triggers. The filtered output 
5 restores the approximate shape and size of the input pulses. The amount of averaging 

increases as the motion artifact increases as determined by the Ensemble _Averaging_Weight 
received from Signal Metrics subsystem. 

[0074] Figure 4 is an example of the Ensemble Averaging subsystem's variable weighting. 
R-Wave triggers mark the start of each pulse. The IR Input line shows pulses corrupted by 
10 motion artifact. The Filtered IR line is a composite of multiple pulses that largely restores the 
original pulse size and shape. Filtered IR is delayed by one second from the IR input. 

[0075] Lowpass Filter - The IR and Red input waveforms are IIR filtered using a weight 
{LPF Weight) received from the Signal Metrics subsystem as follows: 

Lowpass _Waveforms t = Lowpass _Waveforms t -\ + LPF_Weight * {Input ^Waveforms t + 
1 5 Lowpass _ Waveforms h \) 

[0076] During subsystem initialization, the weight is set to a default of 1.0. 

[0077] One-Second Delay Buffer - The IR and Red Lowpass _Waveforms (along with their 
associated age and status) and the RWaveJDccurred input are stored in one-second long 
buffers {IR Inputs, Red_Inputs y Input Valid, Age_Inputs t RWave Inputs). 

20 [0078] Timestamp - The timestamp {Current Time _ctr) is simply a 32-bit counter that is 
initialized to 0 and then incremented every sample period. Several received values are stored 
along with their timestamps in order to reconstruct their age. 

[0079] Ensemble Weight Buffer - The last four Ensemble _Averaging_Weights received 
and their timestamps are stored in the Ensemble _Aver age Buffer and the 
25 Ensemble _WeightJTimestamp. This allows every entry in the one-second delay buffers to be 
associated with its correct weight. Whenever delayed inputs are retrieved from the one- 
second delay buffers, the associated ensemble weight is set to the oldest weight in this buffer 
that is at least as recent as the delayed sample (weight timestamp + one second > 
Cur r en t_ Tim ectr) . 
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[0080] Trigger Qualification - The subsystem selects its trigger from one of two inputs: 
RWavejOccurred or Optical Period. RWavejOccurred is the default trigger. Before 
RWave Occurred is used as a trigger, it passes through the one-second-delay buffer 
mentioned earlier in order to be synchronous with the delayed ER and red samples. Triggers 
5 based on Optical Period input from the Pulse Rate Calculation subsystem are only qualified 
when R-Wave triggers have not been received for at least five seconds. This waiting period 
is deemed sufficient to determine that R-Wave triggers are unavailable. Then, the first trigger 
from the Optical Period is delayed until the first zero-crossing of the subsystem's 
Optical Trigger _WF input waveform. Subsequently, the trigger is derived solely from 
1 0 Optical Period. 

[0081] Ensemble Averaging Model - The subsystem receives IR and Red input samples 
and a trigger to indicate the start of each pulse. Note that the IR and Red inputs have already 
been lowpass filtered and passed through the one-second-delay buffer. It estimates the 
pulsatile component of its current input sample by averaging the ith sample of the current 
15 pulse with the zth sample of the previous estimated pulse. 

[0082] The filter output is therefore a composite of multiple pulses, calculated as each 
sample is received. The filter has the frequency response of a comb filter that only passes 
frequencies at or near the pulse rate and its harmonics. The amount of averaging determines 
the width of the comb filter's "teeth". 

20 [0083] The filter assumes that the zth sample value of the current pulse is roughly equal to 
the /th sample of the previous pulse. Note that i is a ramp function of time, t, that has a 
value of zero at each qualified trigger and increments at each subsequent sample. 

[0084] Ensemble Averaging Filter Equations and Intermediate Variables - The 

Ensemble Averaging filter uses the Ensemble _Aver 'aging Weight, w, which is supplied by the 
25 Signal Metrics subsystem. The following equation shows the basic steps that must be 
performed at the ith sample of every pulse: 

z, = z; + w,(* f -z;) 

where zj denotes the value of z, from one pulse ago. z i is the ensemble-filtered output of the 

subsystem. All z/s are stored in the Composite Pulse Buffer. x t is the most recent output 
30 sample from the one-second-delay buffer. All the variables in equation 1 are scalars. 
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[0085] Composite Pulse Buffers - The composite IR and Red pulses are stored in separate 
Composite Pulse Buffers. The index, /, is reset to the start of the buffers when a trigger is 
received and then incremented for each sample. 

[0086] The buffers must be long enough to store one 20 BPM pulse plus a 10 percent 
5 margin to allow for pulse-rate variability. Therefore, the buffers are updated with at least 1.1 
consecutive composite pulses at 20 BPM and two composite pulses at most pulse rates. If i 
should go past the end of the buffer, normal processing must be suspended until the next 
trigger, during which time the filter's outputs will be set identical to its inputs. 

[0087] When the interval between triggers gets longer, the Composite Pulse Buffers may 
10 not contain recent samples to average at the end of a pulse. 

[0088] The jth sample of the composite pulse is calculated each sample, where j=i+m and 
m is the number of samples between the current and previous qualified triggers, i.e the pulse 
period. Note that j, like /, is a ramp function of t. 

[0089] For j, equation (1) is modified to read: 

Zg - z'g + w t (x t -z') 

15 3 3 ' ' 3 (lb) 

[0090] The same value of w is used to calculate z, and zy. When j reaches the end of the 
buffer, calculation of z, must be suspended until j becomes valid again. 

[0091] Changing pulse periods might cause small discontinuities in the subsystem output 
after each trigger is received. Special processing during the first four samples after a trigger 
20 reduces this effect. The filtered output sample is interpolated between the standard output 
(equation 1) and the second composite pulse (equation lb), so that the filtered output is 
respectively 80%, 60%, 40%, and 20% of z y on these samples. If the second composite pulse 
is not available, the input waveform is used in its place. 

[0092] Initialization, Reinitialization, Clearing and Ignored Samples - Until two 
25 triggers plus one second samples have been received,y is meaningless and z 7 is not calculated. 

[0093] The subsystem contains two methods to recover from an interruption in processing. 
When more than five seconds elapse between triggers, the subsystem is "cleared" at the next 
qualified trigger. The clearing operation sets all of the subsystem's persistent variables to 
their initial values, with the exception of the following variables that are essential to trigger 
30 qualification and maintaining the one-second-delay buffers: 
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1 . The subsystem's current timestamp. 

2. The elapsed times since the last R-wave trigger and the last qualified trigger. 

3. The previous Optical _Trigger_Waveform sample, used to detect zero crossings. 

4. The state of the Trigger Qualification State Machine. 
5 5. One-second delay input buffers 

6. Ensemble Aver aging_ Weight buffer 

[0094] This "clearing" operation is performed because this condition should be rare and 
may indicate a lengthy interruption in processing. The filter outputs are also identical to the 
filter inputs whenever i overflows the end of the Composite Pulse buffer. 

1 0 [0095] When a sample is ignored (not processed) due to lack of valid input for periods 
under five seconds, the Composite Pulse Buffer is left unchanged and the subsystem's 
outputs are marked Invalid. When processing resumes and if Optical ^Period is non-zero, the 
buffer's indices are reset to what they should have been an integral number of pulse periods 
ago and the one-second-delay buffers are reset. This is done because the Composite Pulse 

15 Buffer should still contain an accurate representation of current pulses after a brief 

interruption in processing. The subsystem is "cleared" if the Optical _Period estimate is zero 
(invalid) or samples are ignored for at least five seconds. Most "ignored samples" are 
expected to be due to the adjustments in the oximeter's LED brightness or amplifier gain, 
which generally take less than two seconds. The equations for resetting the indices are as 

20 follows: 

Elapsed ^Periods = m\{Interruption_Duration I Optical JPeriod) 

Samples In Fr actional JPeriod = Interruption Duration - Elapsed Periods * 
Optical Period 

New_Idx = round(/_/dx + Samples In _Fr actional Period) or if Newjdx > Optical ^Period, 

25 New_Idx = round(/ _Idx + Samples Jn_Fr actional _Period - Optical JPeriod) 

j = New_Idx +y - i_Idx, provided j has not already overflowed, and would not overflow, the Composite 
Pulse Buffer. 

ijdx = Newjdx, provided ijdx has not already overflowed, and would not overflow, the 
Composite Pulse Buffer. 
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[0096] The Ensemble Averaging subsystem is reinitialized when the Pulse Lost and Sensor 
Off Detection subsystem determines that a pulse has been re-acquired after being absent for a 
prolonged period of time, or when a sensor is connected. "Reinitialization" means that ALL 
of the subsystem's persistent variables are set to their initial values. This is done because 
5 both of the events that invoke reinitialization make it likely that the Ensemble Averaging 
subsystem's previous pulse representation is no longer current. 

[0097] The elapsed time intervals since the last R-wave and the last qualified trigger are 
incremented on ignored samples. R- Waves or zero-crossings that occur during ignored 
samples are not used to qualify triggers. 

10 [0098] The state transition diagrams of Figs. 5-7 show the state machines for updating each 
of the two pulses maintained in the Composite Pulse Buffers and their indices. 

[0099] Age Metric - The subsystem receives the age in samples, Agejnputs h of its IR and 
Red inputs and outputs the age, AgeJDut^ of its composite pulse outputs. For each sample, 
AgeJDut i, is incremented by the elapsed time since it was last updated and then averaged 
15 with Agejnputs h using the filter weight (Ensemble_AveragingJVeighi), w, used to update 
the composite waveforms. The formula for Age_Out j is: 

age _ out { = age _ out] + m + w t ((age _ in t + N) - (age _ out] +m)) 

where m is the number of samples since Age_OuU was last updated and N is the number of 
samples in one second. The timestamp at which each AgeJDut) is updated must be stored in 
20 order to calculate m. 

[0100] The subsystem must also update AgeJDutj. The formula for AgeJDutj is as above, 
except that j is substituted for /. 

[0101] The subsystem must increment a current timestamp (Current _Time_Ctr) every 
sample, including ignored samples. 

25 [0102] When the subsystem is cleared or reinitialized, all entries in the AgeJDut buffer are 
reinitialized to Age lnputs t . Furthermore, all entries in the buffer for the timestamps at which 
the AgeJDut entries were updated are set to the current timestamp minus 1. These steps are 
taken to assure that the AgeJDut values are no older than the elapsed time since the 
subsystem was cleared or reinitialized. 

14 



[0103] The Ensemble Averaging instance whose waveforms are used for pulse 
qualification and pulse rate uses Rate LPF Weight, which depends solely on frequency 
content. The Ensemble Averaging instance whose waveforms are used to calculate ratio-of- 
ratios and saturation uses Sat LP F Weight, which also depends on whether the 
5 RoR_Variance metric would be better (lower) with the addition of lowpass filtering. These 
weights will range between 0.1 and 1.0, and will not increase more than 0.05 in any single 
step. 

ENSEMBLE AVERAGING WEIGHTS 

[0104] When the subsystem is notified that the Pulse Identification and Qualification 
10 subsystem has just completed evaluation of a potential pulse, the subsystem updates 

ensemble-averaging weights, used by the instances of the Ensemble Averaging subsystem. 
Separate weights are computed for the two Ensemble Averaging instances whose outputs are 
used in computing saturation and pulse rate. These weights are based in part on metrics 
provided by the instance of the Pulse Identification and Qualification subsystem whose input 
1 5 waveforms have NOT been ensemble averaged. 

[0105] The equations for Sat Ensemble _Averaging Filter Weight are as follows: 
x = m&x(Short_RoR Variance, Pulse _Qual_RoR_Variance I 1.5) * 
max(Long_Term_Pulse_Amp_Ratio, 1.0) 
RoR_Variance_Based_Filt_Wt = 0.5 * 0.05 / max(0.05, x) 
20 ArrJProb = (Period_Var - 0. 1 * Short _RoR_Variance - 0.09) / (0.25 - 0.09); 
Arr_Min_Filt_Wt_For_Sat = 0.05 + 0.5 * bound(Arr_Prob, 0, 1.0) 
SatJEnsemble_Averaging_Filter_Weight = max(RoR_Variance_Based_Filt_Wt, 
ArrJdin_Filt_Wt_For_Sai) * 

(1.0 + Pulse_Qual_Score) 
25 SatJLnsemble_Averaging_Filter_ Weight = mm(Sat_Ensemble_Averaging_Filter_Weight, 
1.0) 

where bound(a, b, c) denotes min(max(a, b), c). 

[0106] The above equations result in a default weight of 0.5 for low values of the Ratio-of- 
Ratios variances. Short _RoR_Variance and Pulse jQual_RoR_Variance are both computed 
30 over a three-second interval. The interval for Pulse_Qual_RoR_Variance ends with the 
qualification or rejection of the most recent pulse, which would usually include the most 
recent samples. The weight is reduced by high Ratio-of-Ratios variances, and by high values 
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of Long _Term Pulse _Amp_Ratio that would typically indicate motion artifact. 
Arr_Min_Filt_Wt_For_Sat imposes a minimum value on the ensemble-averaging weight 
(range 0.05-0.55) based primarily on PeriodVar, which quantifies the degree of arrhythmia. 
This is done because ensemble-averaging is not effective for pulses having dissimilar periods. 
5 If the most recent pulse received a good Pulse jQual Score, this can increase the maximum 
value of Sat_Ensemble_Averaging_Filter_Weight from 0.5 to 1.0. 

[0107] The equations for Rate ^Ensemble Aver aging_Filter_Weight are as follows: 
Arr_Prob = (Period_Var - 0.07) / (0.20 - 0.07) 
Arr_Min_FiltJVt_For_Rate = 0.05 + 0.5 * bound(Arr_Prob, 0, 1.0) 
10 x = max(RoR_ Variance _Based_Filt_ Wt 9 ArrJAin_Filt_ Wt_For_Rate) * ( 1 .0 + 
Pulse JQual Score) 

if Short _Term Pulse _Amp_Ratio * Long _Term_Pulse Amp Ratio < 1.0 

x = x I Short _Term Pulse Amp _Ratio 
if Avg_Period > 0 

15 x = x * bound(Pulse__Qual_Score * Qualified _PulseJPeriod I AvgPeriod, 1 .0, 3.0) 

Rate_Ensemble_Averaging_Filter_Weight = min(;c, 1.0) 

[0108] These equations differ from the ones for SatJEnsemble_Averaging_Filter_Weight 
as follows: 

a) The thresholds used to compute Arr_Prob are somewhat lower, because it is desirable 
20 that arrhythmic pulses not be obscured by ensemble averaging prior to pulse qualification. 

b) Small values of Short Term Pulse Amp _Ratio typically indicate that motion artifact has 
just subsided, which means that the ensemble-averaging weight may be quickly 
increased. This has been found empirically to be beneficial for pulse qualification, but 
not for ratio-of-ratios filtering and saturation computation. 

25 c) If the heart skips a beat, with or without prior arrhythmia, the longer-than- average 

Qualified Pulse Period that results will increase the ensemble-averaging weight, so as 
not to obscure the skipped beat from subsequent pulse qualification. 

Definitions : 

Data Inputs 

30 [0109] Avg_Period - Average pulse period reported by Pulse Rate Calculation subsystem 
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[0110] Long JTerm Pulse _Amp _Ratio - Quantifies last pulse amplitude compared to 
historic pulse amplitude. Provided by the Pulse Identification and Qualification subsystem. 
Values substantially larger than 1.0 are typically indicative of motion artifact, and result in 
lower Ensemble _Averaging_Filter_ Weights. 

5 [0111] Period_Var - Period-variability metric from the Pulse Identification and 
Qualification subsystem. Used to gauge the extent of arrhythmia. 

[0112] Pulse_QualJloR_Variance - RoR_Variance metric from the Pulse Identification and 
Qualification subsystem. For instance, a value of 0.10 would indicate that the average 
difference between consecutive pulse periods is 10% of Avg Period. 

1 0 [0113] Pulse jQual_Score - Score computed by the pulse qualification neural net in the 

Pulse Identification and Qualification subsystem. Zero is extremely poor and 1.0 is excellent. 

[01 14] Qualified _Pulse_Period - Most recent pulse period qualified by the Pulse 
Identification and Qualification subsystem. 

[0115] Short _Term_Pulse_Amp _Ratio - Quantifies last pulse amplitude compared to 
1 5 previous pulse amplitude. 

Outputs 

[0116] Frequency Ratio - Ratio of Mean_IR_Frequency_Content to pulse rate. 

[0117] LPF RoR Variance - Quantifies variability of ratio-of-ratios. Computed over a 9- 
second window from LPF Scaled _Waveforms. 

20 [0118] Rate_LPF_Weight - Lowpass filter weight to be used by the instance of the 

Ensemble Averaging subsystem that preprocesses waveforms used for pulse qualification and 
pulse rate calculation. 

[0119] RoR_Variance - Quantifies variability of ratio-of-ratios. Computed over a 9-second 
window from Scaled ^Waveforms. A value of 0.10 would indicate that sample-to-sample 
25 ratio-of-ratios values differ from the mean ratio-of-ratios value by an average of 10% of the 
mean ratio-of-ratios value. 

[0120] Sat _Ensemble_Aver aging Filter ^Weight - Ensemble-averaging weight to be used by 
the instance of the Ensemble Averaging subsystem that preprocesses waveforms used for 
pulse qualification and pulse rate calculation. 
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[0121] Sat JLPF Weight - Lowpass filter weight to be used by the instance of the Ensemble 
Averaging subsystem that preprocesses waveforms used for pulse qualification and pulse rate 
calculation. 

[0122] Scaled ^Waveforms - Scaled versions of IR and Red Pre_ProcessedJWaveforms. 

5 [0123] Short JRoR_Variance - Quantifies variability of ratio-of-ratios. Computed over a 3- 
second window from Scaled ^Waveforms. 

Internal Variables 

[0124] Arr_Prob - Likelihood of arrhythmia that would limit the amount of ensemble 
averaging. Based on Period_Var, with threshold that are specific to each of the two 
1 0 Ensemble _Aver aging JFilter ^Weights. 

[0125] Arr_Min_Filt_Wt_For_Rate, Arr_Min_Filt_Wt_For_Sat - Minimum values for the 
two Ensemble _Averaging_Filter_Weights, based on their respective Arr Prob values. 

[0126] LPF _Scaled_Waveforms - Lowpass-filtered version of Scaled ^Waveforms, used to 
compute LPFRoRVariance. 

1 5 [0127] Mean_IR_Frequency Content - Estimate of the mean frequency content of the IR 
input waveform. 

[0128] RoR_ Variance Based _Filt_ Wt - Component for 

Ensemble _Averaging_Filter _Weights based on RoR Variance metrics and 

Long_Term_Pulse_Amp_Ratio. 
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